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Topological magnetoelectric effect in a three-dimensional topological insulator is a novel phe¬ 
nomenon, where an electric field induces a magnetic field in the same direction, with a universal 
coefficient of proportionality quantized in units of e 2 /2 h. Here we propose that the topological mag¬ 
netoelectric effect can be realized in the zero-plateau quantum anomalous Hall state of magnetic 
topological insulators or ferromagnet-topological insulator heterostructure. The finite-size effect is 
also studied numerically, where the magnetoelectric coefficient is shown to converge to a quantized 
value when the thickness of topological insulator film increases. We further propose a device setup 
to eliminate the non-topological contributions from the side surface. 

PACS numbers: 73.43.-f 73.20.-r 85.75.-d 


The search for topological quantum phenomena has 
become an important goal in condensed matter physics. 
Topological phenomena in the physical systems are de¬ 
termined by topological structures and are thus univer¬ 
sal and robust against perturbations, and the electro¬ 
magnetic response is usually exactly quantized [1]. Two 
well-known examples of topological quantum phenom¬ 
ena are the flux quantization in superconductors [2] and 
Hall conductance quantization in the quantum Hall effect 
(QHE) [3]. The remarkable observation of such topolog¬ 
ical phenomena is that the quantization is exact, which 
provide the precise values of fundamental physics con¬ 
stants such as Plank’s constant h [4]. 

The recent discovery of the time-reversal (T) invariant 
(TRI) topological insulator (TI) brings the opportunity 
to study a large family of new topological phenomena [5, 
6]. The electromagnetic response of a three-dimensional 
(3D) insulator is described by the topological 6 term [7-9] 
of the form 

Op 2 !' 

Sg = —- d 3x( itE-B, (1) 

together with the ordinary Maxwell terms. Here E and B 
are the conventional electromagnetic fields inside the in¬ 
sulator, e is the charge of an electron, and 0 is the dimen¬ 
sionless pseudoscalar parameter describing the insulator, 
which refers to the axion field in particle physics [10]. Un¬ 
der the periodic boundary condition, all physical quan¬ 
tities are invariant if 6 is shifted by integer multiples of 
27 t. Therefore, all TRI insulators are described by either 
6 = 0 or 0 = n (modulo 27r). TIs are defined by 0 = 7r, 
which cannot be connected continuously to trivial insu¬ 
lators, defined by 0 = 0, by TRI perturbations. With 
open boundary condition, the effective action is reduced 
to a (2+l)D Chern-Simons term on the surface, which 
describes a surface QHE with half-quantized surface Hall 
conductance [7]. Such a topological 6 term with a uni¬ 
versal value of 0 = n in TIs leads to a magnetoelectric 


effect with coefficient quantized in units of e 2 /2 h, known 
as the topological magnetoelectric effect (TME), i.e., an 
electric field can induce a magnetic polarization, whereas 
a magnetic field can induce an electric polarization. To 
obtain the quantized TME in TIs, as is first suggested in 
Ref. [7], one must fulfil the following stringent require¬ 
ments. First, introduce a T-breaking surface gap by fer¬ 
romagnetic (FM) ordering, where the magnetization of 
FM points inward or outward from the surface. Second, 
finely tune the Fermi level into the magnetically induced 
surface gap and keep the bulk truly insulating. Third, the 
film of TI material should be thick enough to eliminate 
the finite-size effect, therefore the TME is exactly quan¬ 
tized. Several other theoretical proposals [11-14] have 
been made to realize the TME; however, observing the 
TME in TIs experimentally is still challenging. 

In this paper, we propose to realize TME effect in 
the newly discovered quantum anomalous Hall (QAH) 
state [15, 16]. Recently, a new zero-plateau QAH state 
in a magnetic TI has been theoretically predicted [17] 
and experimentally realized [18, 19]. The magnetic TI 
studied in the QAH experiment develops robust FM at 
low temperature. In the magnetized states, the magnetic 
domains of the material are aligned to the same direction, 
and the system is in a QAH state with a single chiral edge 
state propagating along the sample boundary, where the 
Hall conductance a xy is quantized to be ±e 2 //i. The 
zero-plateau state, on the contrary, appears around the 
coercivity when the magnetic domains reverse, where a xy 
shows a well-defined zero-plateau over a range of mag¬ 
netic field around coercivity while longitudinal conduc¬ 
tance a xx —» 0 [shown in Fig. la]. In such a state, the 
Fermi level is in the magnetization induced surface gap, 
fulfilling the first two conditions above, providing a good 
platform to observe the TME effect as we will discuss in 
details below. However, due to finite thickness in mag¬ 
netic TI, the TME is non-quantized. Therefore, we fur- 
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FIG. 1. (color online). Zero-plateau QAH state and magnetic 
field dependence of a xy . (a) Magnetic TI in an external field 
H, and sketch of o xy as a function of H. a xy = 0 plateau oc¬ 
curs at the coercivity. The arrow indicates the magnetization 
direction, (b) FM-TI-FM heterostructure, and a xy vs H. The 
zero-plateau appears in the hysteresis loop due to different iff 
and H 2 of FM A and B. The blue and red arrows indicate the 
magnetization direction of FM A and B, respectively. 


ther propose to realize the quantized TME effect in the 
zero-plateau QAH state of the FM-TI-FM heterostruc¬ 
ture as shown in Fig. lb, where an in-plane ac magnetic 
field induces an electric current in the same direction, or 
an in-plane electric field induces a magnetic field. The 
finite-size effect is also studied numerically, where the 
TME coefficient is shown to converge to a quantized value 
when the thickness of TI film increases. Finally, we pro¬ 
pose a device setup where the non-topological contribu¬ 
tion from the side surface is negligible. 

The TME described by the topological 6 term implies 
that a quantized magnetic polarization is induced by an 
electric field, given by 

0 e 2 

m = < 2 > 

Such response can be understood in terms of a surface 
Dirac fermion picture. With antiparallel magnetization 
of the two FM layers as shown Fig. 2b, an in-plane electric 
field E = E y y induces the Hall currents J* = a* y z x E 
and J fa = cr xy z x E on top (z = d/2 and denoted as 
superscript t) and bottom (b and z = —d/2) surfaces, 
respectively. Since the surface massive Dirac fermion 
gives rise to half-integer Hall conductance a xy = — a b y = 
(9/2tt)( e 2 /h), the currents J* = —J b are opposite and 
form a circulating total current. Inside the sample, such 
a circulating current can be viewed equivalently as the 
surface bound current generated by a constant magne¬ 
tization M = — J f y = — (9/2n)(e 2 /h)'E. Therefore, the 
TME essentially originates from the half-quantized sur¬ 
face Hall conductance. 

First we examine theoretically the TME in the zero- 
plateau QAH state observed in experiments. The low- 


energy physics of this system consists of Dirac-type sur¬ 
face states only [17, 20]. At the coercivity, both ran¬ 
dom magnetic domains that formed in the sample, and 
the exchange field A introduced by the FM ordering are 
spatially inhomogeneous. The 3D spatial average of ex¬ 
change field vanishes ((A) ave = 0). However, due to 
an unavoidable top-bottom asymmetry, the top and bot¬ 
tom surfaces may feel an opposite nonzero exchange field 
A* = —A & = Aq. In this case, the zero-plateau state 
is described by the mean field effective model which has 
the generic form as 7f 0 = k y a\ ® T 3 — k x 02 <S> 73 + A 0 O 3 G> 
73 + ml <g> T\. with the basis of 1 1 t), 1 1 J,), |b t) and 
\b |), where f / 4 - represent the spin up/down states, 
respectively, cr, and r,; (i = 1,2,3) are Pauli matrices 
acting on spin and layer, respectively, m describes the 
hybridization between the top and bottom surface states. 
If m = 0, due to the opposite half-integer Hall conduc¬ 
tance contributions from the top and bottom surfaces 
a xy = ~ a xy — sgn(A 0 )(e 2 /2 h), the system has a x ° y = 0, 
which gives rise to a quantized TME as discussed above. 
However, a nonzero m will mix the circulating current 
J* and J 6 , therefore, the TME is no longer quantized. 
In reality, the exchange field depends very much on the 
microscopic details of the randomness in magnetic do¬ 
mains. However, we emphasize that the TME of the 
zero-plateau state in a magnetic TI is in general nonzero 
and non-quantized. 

To realize the quantized TME effect, the TI film should 
be thick enough so that the hybridization between top 
and bottom surfaces is negligible. Therefore, we pro¬ 
pose that a quantized TME can be realized in the zero- 
plateau state of FM-TI-FM structure as shown in Fig. lb. 
The FM insulators A and B have different coercivity Hf 
and if I, respectively. Assume that both FM A and B 
have an out-of-plane magnetic easy axis, and the same 
sign of the exchange coupling parameter to TI surface 
states. When A and B have antiparallel magnetization, 
the system is in a zero-plateau QAH state with o xy = 0, 
which is contributed by a% y + a b y as (1/2 — 1/2 )(e 2 /h) 
or (—1/2 + 1/2 )(e 2 /h). Such a magnetization configura¬ 
tion can be easily achieved in the hysteresis loop by an 
external field if with if/ < |if| < if/, and then remove 
if. Experimentally, to achieve the TME in this setup, 
a good proximity between FM and TI is necessary. TI 
material can be chosen as Bi y Sbi_ y Te 3 , where the Dirac 
cone of the surface states is observed to be located in 
the bulk band gap [21]. The candidate FM materials 
are Cr 2 Ge 2 Te 6 (CGT), Cr^BqSb^-xTes (CBST) with 
0.3 < x < 0.46 and V a ,(Bi,Sb) 2 - a ,Te 3 (VBST). All of 
them are FM insulators with an out-of-plane easy axis, 
and have good lattice match with Bi 2 Te 3 family mate¬ 
rials. CGT is a soft FM insulator with T c ~ 61 K and 
if c < 100 Oe [22], and it also shows good proximity with 
Bi 2 Te 3 [23]. CBST with 0.3 < x < 0.46 is a FM insulator 
with T c = 40-90 K and H c ~ 1.0 x 10 3 Oe [24], VBST 
with 0.1 < x < 0.3 is a FM insulator with T c = 30-100 K 
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FIG. 2. (color online). TME effect, (a) Illustration of the 
ac electric current induced by an ac magnetic field B x , the 
current density jfP{z) is defined in Eq. (5). (b) Magnetic 
field By induced by applying an electric field E y through a 
capacitor. E y (with direction into the paper) will induce Hall 
currents J‘ and J 6 for antiparallel magnetization, (c) The 
functions k(z) and r/(z) for different thickness 6, 10, 20 QL. 
Each QL is about 1 nm thick. 


and H c ~ 1.0 x 10 4 Oe [25], 

TME As we discussed previously, an electric field will 
induce a topological contribution to bulk magnetization. 
From the constituent equation H = B /n — M, with 
H = 0 and B continuous, we have on the middle of 
side surface (parallel to z) B = —fi(e 2 /2h)E y y. Here 
H is the material-dependent magnetic permeability. Tak¬ 
ing fi ss hq, E y = 10 5 V/m, we get the magnitude of 
magnetic field 2.43 x 10 ~ 6 T, which is easily detectable 
by present superconducting quantum interference devices 
(SQUID). The stray magnetic field effect can be well sep¬ 
arated from the quantized TME by ac modulation of the 
electric field and phase-locking detection, where the ac 
frequency is quasi-static around 10-100 Hz. Moreover, a 
gradiometer sensor in SQUID could also screen the ho¬ 
mogeneous stray field. 

The TME also indicates the induction of a parallel 
polarization current when an ac magnetic held is ap¬ 
plied. Consider the process of applying an ac mag¬ 
netic held B = B x x. as shown in Fig. 2a. A circulat¬ 
ing electric held E parallel to side surface (parallel to 
B) is generated due to Faraday’s law of induction, where 
E 4 = —E b = (dB x /dt)(d/ 2)y. Such an electric held will 
induce a Hall current density j 2D = j 2D + j 2D , where 
j 2D = a xy z x E f and j 2D = a xy z x E b . Therefore the 
total current J = j 2D £ = Jx, where 


J 


9 e 2 dB x 
7 t 2 h dt 


Id. 


( 3 ) 


Here d and t are the thickness and width of the TI him as 
shown in Fig. 2a, and 9 —> n when d is large enough. For 
an estimation, take B x = B 0 e~ lult ^ B 0 = 10 G, lo/2tt = 
1 GHz, d = 20 nm, 9/ir ss 0.91 (hnite-size effect taken 
into account as in Fig. 3), and i = 500 /jm, we have J = 



FIG. 3. (color online). Finite-size effect of TME. (a) The 7 
and k(z = 0) as a function of d. The inset shows 7 plotted 
vs the inverse of thickness 1 /d. (b) The current amplitude Jo 
scales linearly with d. Here 6/n = 7 . 


—iJoe~ lult with J 0 = 1.11 nA, in the range accessible by 
transport experiments. Moreover, as shown in Fig. 3b, 
the current amplitude Jo scales linearly with thickness 
d , for 9 is a linear function of 1/d with (1 — 9/n) oc 1/d, 
i.e., thicker him gives rise to larger TME. 

Finite-size effect Due to the hnite-size conhnement 
along z direction, the TME effect is not quantized when 
the TI him is thin. However, as we shall show be¬ 
low, the TME coefficient converges quickly into the 
quantized value as the him thickness d increases. The 
generic Hamiltonian of a TI thin him can be written as 
TL 2 D(k) = f^ 2 /2 dz'H 3D (k, z). Here k = (k x ,k y ), and 
we impose periodic boundary conditions in both x and y 
directions. The magnetoelectric response of such a thin 
him can be directly calculated with the Kubo formula. 
With the 3D in-plane current density operator defined as 
j 3 D (k, z) = (e/ 7 i) 9 k% 3 D(k, z), we can write down a dc 
current correlation function 


n xy (z,z') 


x 21 m 


= / d2k E /('»*> 

n^m 

(^nk| jx (k; z) | Um]f) (u m k| j y (k, Z ) |u n k) 
(Cnk Cmk)“ 


where \u n \f) is the normalized Bloch wavefunction in 
the n-th electron subband satisfying TL 2 D(k) |u„k) = 
enk|u„k) , and /(e) is the Fermi-Dirac distribution func¬ 
tion. The Kubo formula for magnetic held B y induced 
by a uniform external electric held E y is then given by 


B y (z) = -^k(z)E v , (4) 

where k(z) = dz x sgn(z - zi) J d/ f /2 dz 2 B xy {z 1 ,z 2 ) 
is a dimensionless function. Here sgn(z) gives the sign 
of z. Similarly, the current density j 3D induced by a 
uniform external ac magnetic held B x of frequency w/27t 
is given by 


E (z) = -iu>—r)(z)B ; 


( 5 ) 
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where 77(2) = 2 J^ 2 dzi ZiH xy (z, Zi) is also dimension¬ 
less. 

The response formulas above are generic for any TI 
system and do not rely on a specific model. For con¬ 
creteness, we adopt the effective Hamiltonian in Ref. [ 26 ] 
to describe the low-energy bands of Bi2Te3 family mate¬ 
rials, "H3D(k, z) = el®l + d 1 T\®l-\-d 2 T2®(J3 + d ? ’T3®l- 
A (z)t 3 (g) cr 3 + iAid z r 2 g> <J 2 . Here Tj and aj (j = 1 , 2 , 3 ) 
are Pauli matrices, e(k, z) = — Did 2 + D 2 {k 2 + k 2 ), 
d 1 ' 2 ’ 3 (k^) = C A 2 k x ,A 2 k y ,B 0 - B x d 2 z + B 2 (k 2 x + k 2 y )), 
and A(2) is the 2-dependent exchange field. We then 
discretize it into a tight-binding model along 2-axis be¬ 
tween neighboring quintuple layers (QL) from TLsDi and 
assume A (z) takes the values ±A S in the top and bottom 
layers, respectively, and zero elsewhere. Fig. 2 c shows the 
numerical calculations of k(z) and 77(2) for thin films of 
6, 10 and 20 QL, where we set a typical surface exchange 
field A s = 50 meV. All the other parameters are taken 
from Ref. [ 27 ] for (Bio.iSbo.g)2Te 3 . The bulk value of 
k(z) at 2 = 0 as a function of d is plotted as black line in 
Fig. 3 . As is consistent with the topological field theory, 
k(z) in the bulk tends to 1 and becomes quantized as the 
thickness d increases, whereas 77(2) is bounded within a 
finite penetration depth to the top and bottom surfaces. 
The shape of functions k(z) and 77(2) near surfaces re¬ 
main almost unchanged as the thickness d varies. 

To characterize the deviation from topological quan¬ 
tization of TME in TI thin films, we further define the 
dimensionless number 

\ rd/2 1 rd./2 

7 = - / dz k(z) = - dz 77 ( 2 ) , ( 6 ) 

“ J-d/2 « J — d/2 

which is the mean value of k(z) or 77 ( 2 ) (which are 
equal to each other). The average magnetic field in re¬ 
sponse to the external electric field E y is then R" lean = 
— 7 /i(e 2 / 2 h)E y , whereas the total 2D current density in¬ 
duced by external magnetic field B x is given by j 2D = 
f dzj^ D (z) = — iuo r )d{e 2 /2h)B x . Compared to Eq. (3), 
we get 9/tt = 7 . The value of 7 as a function of d 
is shown as the red line in Fig. 3, where 7 —>■ 1 with 
d —> 00 . This shows the TME effect is quantized as 
the system is in the thermodynamic limit. In fact, as 
shown in the inset of Fig. 3, the value of 1 — 7 scales lin¬ 
early with 1/d as the thickness d —> 00 , which indicates 
•f—d /2 — K ( z )) = const, when d is large enough. This 

is simply because the function 1 — k(z ) is nonzero only 
near the top and bottom surfaces, and its shape is inde¬ 
pendent of the thickness d , as shown in Fig. 2c. Since 
(1 — 7 ) oc 1/d, Jo (xdd = 77 r d is a linear function of d, 
as is shown in Fig. 3b. 

Discussion The TME effect in the setup Fig. 2 is not 
quantized when the side surface (parallel to z) is not 
gapped. The gapless side surface states may give rise to 
a non-topological contributions to the TME effect. To 
eliminate such non-topological contribution, one can use 
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FIG. 4. (color online). Schematic of FM-TI-FM heterostruc¬ 
ture to observe the quantized TME, where the side surface of 
TI is gapped by FM proximity in (a). Au is the electrode. 
Such geometry of TI can be made by lithography. In reality, 
such configuration cook in experiments may like (b), where 
the side surface is uneven and step-like. In this case, the side 
surface is gapped by both FM ordering and quantum confine¬ 
ment. 


the device setup as shown in Fig. 4, where the side surface 
is gapped either by FM order or quantum confinement. 
Also in this setup, if FM A and B have opposite sign of 
exchange coupling parameter to TI surface, only paral¬ 
lel magnetization is needed to realize TME. Recently, the 
surface QHE has been realized in bulk Bi 2 _ a: Sb a ,Te 3 _ySej / 
TI [28, 29], where the systems exhibit surface-dominated 
conduction even at temperatures close to the room tem¬ 
perature, whereas the bulk conduction is negligible. Such 
experimental progress on the material growth and rich 
material choice of TI and FM insulator make the realiza¬ 
tion of the quantized TME in TIs feasible. 
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Note added: During the preparation of our manuscript, 
we learned of an independent work on a similar prob¬ 
lem [30]. However, their experimental proposal to ob¬ 
serve the TME is different from our results. 
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